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ABSTRACT 



Aims. We aim at resolving the circumstellar environment around /3 Pic in the near-infrared in order to study the inner 
planetary system (< 200 mas, i.e., ~4 AU). 

Methods. Precise interferometric fringe visibility measurements were obtained over seven spectral channels dispersed 
across the H band with the four-telescope VLTI/PIONIER interferometer. Thorough analysis of interferometric data 
was performed to measure the stellar angular diameter and to search for circumstellar material. 

Results. We detected near-infrared circumstellar emission around /3 Pic that accounts for 1.37% ± 0.16% of the near- 
infrared stellar flux and that is located within the field-of-view of PIONIER (i.e., ~200mas in radius). The flux ratio 
between this excess and the photosphere emission is shown to be stable over a period of 1 year and to vary only weakly 
across the H band, suggesting that the source is either very hot (> 1500 K) or dominated by the scattering of the stellar 
flux. In addition, we derive the limb-darkened angular diameter of /? Pic with an unprecedented accuracy (6'ld= 0.736 
± 0.019 mas). 

Conclusions. The presence of a small H-band excess originating in the vicinity of /3 Pic is revealed for the first time 
thanks to the high- precision visibilities enabled by VLTI/PIONIER. This excess emission is likely due to the scattering 
of stellar light by circumstellar dust and/or the thermal emission from a yet unknown population of hot dust, although 
hot gas emitting in the continuum cannot be firmly excluded. 

Key words. Instrumentation: high angular resolution — techniques: interferometric — circumstellar matter 



1. Introduction 



The young (~12t| Mvr. IZuckerman et al.|[200l[ ) 



A6V-type 



>. star 13 Pic (HD 39060, A6V, 19.3 pc) is surrounded by a fa- AU) 



X 



et al 
spect 



mous planetary system, which is a prime target "or under- 
standing pla netary system formatio n and evolut: on. Since 
.its discovery (jSmith fc Terrilelll984[ ). successive generations 
"of telescopes have reported the detection of aii edge-on 
debris disc with several distinctive features suggestive of 
a multiple-belt system (|Telesco et al.l 20051). st ir-grazing 
comets ("falling eva porating bodies". iBeust et al. 1990), cir- 
cumstellar gas (e.g.. lHobbs et al-lligSa Roberge et'aLl l2006). 
and a 9-Mjiip planetary companion orbiting at a projecte d 
distance of approximately 4.3 AU (iLagraiige et 311120091) 



to the outer disc ( Lagrange et al.ir2012aD . 
Over the past few years, the close environment (< a few 
)f (3 Pic has been the focus of several stud ies trying to 



120101 : 



120011 : iMouillet et aLlll997D inchned by ~4° with 



detec t a putative sub-stellar companion (e.g., lAbsil et al.l 



iLagrange et aLll2012b[ ). In this paper, we present 



The existence of other planets seems likely (|Freisb itter et ah! 
[2007) and might explain several asymmetries id ;ntified in 
the debris disc, including a warp at ~50 AU (jAugereau 



high- iccuracy H-band interferometric observations of /3 Pic 
that were originally aimed at a deep search for faint com- 
panic ns in the inner disc but which revealed instead the 
presence of a more extended circumstellar emission. Us- 
ing our spectrally dispersed measurements and revisiting 
archit^al K-band VLTI/VINCI measurements, we constrain 
the r ature of this excess emission and discuss briefly its 
origii I. 

2. O Dservations and data reduction 



* Based on data collected at the ESO Paranal Observatory 
under commissioning time and programme 088.C-0266. 
** FNRS Postdoctoral Researcher. 



Interferometric data were obtained with the PIONIER com- 
biner (Lc Bouquin et al. 2011) at the VLT interferometer on 
three different occasions (see Table [T]). The fringe measure- 
ments were dispersed over seven spectral channels across 
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Fig. 1. Sampling of tiie Fourier {u,v) plane obtained for 
the complete data set (one colour per spectra l channel). The 
orientation of the outer disc midplane (i.e., 29?5 lBoccaletti et al.l 
l2009h is represented by the black dashed line. 

Table 1. Overview of the data obtained with VLTI/PIONIER 
on /? Pic. 



Fig. 2. Expected squared visibility of the limb darkened pho- 
tosphere (blue solid line) as a function of the spatial frequency, 
along with the measured squared visibilities and related 1 — a 
error bars (one colour per wavelength). The thickness of the 
blue solid line corresponds to the 3-a error related to the uncer- 
tainty on the stellar diameter. The best-fit model is represented 
by the dotted blue line with the residuals of the fit given in the 
middle panel. It corresponds to a limb-darkened photosphere of 
0.736 ± 0.015 ± 0.012 mas in diameter surrounded by a uniform 
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2010/12/20 
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1,2,3 
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2011/11/02 


DO-Gl-HO-Il 


0.9" 


14 


2,3,4 



" An observing block consists in a single observation (sec more details in 
ILe Bouguin et al.ll2011l) . 

Calibrator stars correspond to (1) HD 34642 (KlIV), (2) HD 35765 
(KlIII), (3) HD 39640 (G8III), (4) HD 46365 (K3HI), and (5) HD 223825 
(G9HI). 

the H band (1.50 - 1.80 ^m) but in the following we discard 
the first spectral channel (1.50 - 1.55 /im) that is systemati- 
cally of lower quality and not suitable for the high-precision 
visibility measurements required by our programme. Ob- 
servations of /3 Pic were interleaved with observations of 
reference stars to calibrate the instrumental contribution 
in the observed quantities. Calibrators were chosen close 
to /3Pic, in ter ms of both po s ition and magnitude, from 
the catalogue of iMerand et al.l (|2005D . The total u-v plane 
covered by the observations is shown in Fig.ffl 

Data were reduced and c alibrated with the pndrs pack- 
age (jLe Bouauin et aLll201ll ) . We focus here on the squared 
visibilities (V^) to measure both the stellar angular diam- 
eter and search for circumstellar material. The final cal- 
ibrated data set (V^) is shown in Fig. [5J The search for 
faint companions by means of a closure-phase analysis will 
be presented elsewhere (Lagrange et al. in prep). 



3. Data analysis 

The calibrated were fitted to a range of models consist- 
ing of an oblate limb-darkened photosphere surrounded by 
a uniform emission ("disc") filling the entire field-of-view 
of PIONIER on the auxiliary telescopes. Under typical 
seeing conditions, this field-of-view can be approximated 



by a Gaussi an profile with a full width at half maximum 
of 400 mas (jAbsil et al.ll201lh . equivalent to 4 AU in ra- 
dius at the distance of /3 Pic. Our model is based on two 
free parameters, namely the limb-darkened angular diam- 
eter of the star and the disc/star contrast. The distortion 
of the photo sphere produced b y the rapid rotation (v sini 
= 130km/s. lRover et al.ll2007l ) was considered to produce 
a realisti c model of the sta r. Following the parametric ap- 
proach of lAbsil et al.l ()2008f ). the distorted photosphere was 
modeled by an ellipse with an oblateness of 1.038 and a ro- 
tation axis perpendicular to t he outer disc midplane (which 
has a position angle of 29?5. iBoccaletti et al.|[2009l ). The 
expected from the limb-dar kened photosphere was then 
estimated according to iHanburv Brown et al.l (j 19741 ) con- 
sidering a linear lim b-darkening H-band coefficient of 0.24 
(jClaret et al.lll995r i. 

To search for circumstellar material, we compared the 
measurements to the expected of the stellar photo- 
sphere. The circumstellar emission then appe ars as a 
deficit at short baseline lengths as detailed in lAbsil et al.l 
(120061 ). In this first approach, we assume that the disc/star 
contrast does not depend on the wavelength and fit all mea- 
surements simultaneously with this single parameter. We 
used the bootstrapping method to compute the statistical 
error bars on the stellar angular diameter and the disc/star 
contrast. Considering correlation between the baselines, 
the spectral channels, and the successive measurements of 
the same baseline, we found that correlations mostly oc- 
cur between all spectral channels of the same measurement. 
Hence, we performed the bootstrapping by drawing all spec- 
tral channels together, which corresponds to 126 indepen- 
dent data sets out of our 6 x 6 x 21 = 726 visibility mea- 
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Fig. 3. Normalised map for a model of an oblate pho- 
tosphere (mean limb-darkened angular diameter in horizontal 
axis) surrounded by a uniform disc filling the whole field-of-view 
(in relative flux, vertical axis). The diameter predicted by the 
surface-brightness relations is represented by the vertical lines 
and the position of the best-fit model by a white filled circle. 



surements. We randomly produced 100 data sets from the 
original data and fitted them separately to our model. The 
standard deviation on the derived best-fit parameters then 
gave us reliable error bars. Finally, we considered an ad- 
ditional error to take the chromatic behaviour of the beam 
combiner and the different colours between /3 Pic (A6V) 
and its calibrators (G8III to K4III) into account. Based 
on the spectral shape of the transfer function across the H 
band (maximum variation of 40% over the entire band) , we 
derived a maximum systematic error on the of ±0.2% in 
a single spectral channel (see more details in iDefrere et al.l 
I2OII). This value must be considered as very conservative 
since all other sources of systematic error are expected to be 
much smaller thanks to our optimised observing strategy, 
i.e. spectrally-dispersed observations obtained at various 
epochs on various configurations with different calibrators 
that are close to the science target in both magnitude and 
position. Using the approach described above, we produced 
the normalised map shown in Fig. [3] as a function of the 
meai^B angular diameter of the photosphere and of the H- 
band disc/star contrast. The two contours at 1 and 5-cr in- 
dicate the parameter space that falls within the respective 
confidence levels taking the number of degrees of freedom 
in the ^ distribution into account. In the present case, we 
assume that all spectral channels are fully correlated so that 
there are 21x6 independent measurements and two param- 
eters to fit, hence 124 degre es of freedom. Using the surface 
brightness relations (SBR. iKervella et al.|[20M ) over vari- 
ous wavelength bands (B, V, J, H, K), we derived a stellar 
angular diameter of 0ld — 0-712 ± 0.010 mas. This value 
does not significantly depend on the chosen photometric 
bands, which is a good sign of robustness. For a purely 
photospheric model, the best- fit diameter is 0.875 ± 0.011 
± 0.012 mas, where the second and the third terms repre- 
sent the statistical and the systematic errors, respectively. 
This value is at more than 10-ct from the one given by the 




^ The mean angular diameter corresponds to the geometric 
mean of the minor and major axes of the elliptical photosphere. 
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Fig. 4. Best-fit resolved fluxes given as a function of wave- 
length, assuming that /3 Pic emits as a blackbody of tempera- 
ture 8200 K. The solid lines give the flux expected from black- 
body sources of various equilibrium temperatures and the op- 
timum emitting surface area (computed by minimisation). 
The dashed line represents the flux from a blackbody source of 
8200 K or, equivalently, from scattered stellar light (xr=0.91). 
The K-band flux obtained from archival VLTI/ VINCI measure- 
ments is shown for comparison but is not included in the fit (see 
more information in the main text). 



SBR. Moreover, the residuals show a systematic trend 
versus spatial frequencies (Fig. [51 bottom panel) . Now ac- 
counting for a possible excess emission, the best-fit model is 
obtained for a stellar diameter of 0.736 ± 0.015 ± 0.012 mas 
and a disc/star contrast of 1.37 ± 0.10 ± 0.13%. The best- 
fit diameter is fully compatible with the one given by the 
SBR, and no other trend is seen in the residuals (Fig. 
middle panel). We conclude that our dataset unambigu- 
ously demonstrates the existence of an H-band circumstel- 
lar emission around /3 Pic and in the field-of-view of PIO- 
NIER. This excess emission is repeatedly detected in each 
run with at least a 5-(t confidence level and with compatible 
values within error bars. Using our best-fit angular diame- 
ter, we derived best-fit flux ratios of 1.48 ± 0.13 ± 0.17% 
and 1.32 ± 0.10 ± 0.13% by fitting the data of observing 
runs A&B and C separately. This suggests that the ex- 
cess emission was present and did not vary during the time 
covered by our observations, i.e. approximately one year. 

Finally, we revisited the K-band measurement ob- 
tained in 2002 with the VLTI/ VINCI instrument. These 
observations used a single and relatively long baseline 
length (~92 m), making it impossible to lift the degeneracy 
between the s tellar diameter and a circumstellar emission 
(see Fig. 7 in Idi Folco et all [20041 ). By using our best-fit 
stellar angular diameter (i.e., 0.736 ± 0.015 ± 0.012 mas), 
we derived a best-fit K-band flux ratio of 0.76% ± 0.49%, 
which will be used in the following discussion. 

4. Constraining the resolved emission 

To constrain the origin of the resolved near-infrared emis- 
sion, we first tried to fit the data with d isc models of var- 
ious geometries (see similar analyses in lAbsil et al.ll2009l : 
Defrere et al. 2011). This analysis provided only a con- 
straint on the minimum spatial extent of the near-infrared 
excess emission zone, which must be resolved by PIONIER 
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on all baseline lengths and orientations (i.e., further than 
^0.1 AU from the star). Then, we fitted the data corre- 
sponding to each spectral channel separately, considering 
the procedure described in the previous section and the 
stellar angular diameter derived with all data (i.e., 0.736 ± 
0.015 ± 0.012 mas). Given the low spectral resolution of 
PIONIER (R«40) , we assumed that /3 Pic emits as a black- 
body of temperature 8200 K and converted the dispersed 
best-fit contrasts into fluxes. The best-fit fluxes were then 
fitted to pure blackbody models of various temperatures as 
represented in Fig. 21 This figure shows that a relatively 
high temperature (>1500K) is required to fit the data in 
a satisfactory way {Xr < 5). Including the K-band best- 
fit flux in the fit would reinforce this statement, although 
this must be considered with care given the eight to nine 
years separation and the greater uncertainty of the K-band 
measurements. Models corresponding to temperatures in 
the range of 1500-2000 K would suggest the presence of hot 
material in the field-of-view of PIONIER while higher tem- 
perature models would rather point towards the scattering 
of stellar hght by circumstellar dust (^8200 K, Xr=0-91) or 
towards a binary companion (>2000 K). The latter scenario 
can, however, be excluded with good confidence based on 
other studies (jAbsil et al.ll2010l : iLTgrange et al.ll2012bf ) so 
it is not addressed in the following discussion. 

5. Discussion 

The most attractive scenario for explaining the wavelength- 
independent fiux ratio across the H band is the scattering of 
stellar light in the outer debris disc {Xr = 0.91, see Fig.S]). 
Since the disc is seen edge-on, a significant amount of stel- 
lar light could be forward- (or back-) scattered towards the 
interferometer field-of-view. To assess the validity of this 
scenario, we perform ed radiative tra,nsfer modelling with 
the GraTer software (jAugereau et aLlllQQQT ). This software 
computes synthetic fluxes and images of dust discs in scat- 
tered light and thermal emission over various parameters 
for the dust composition and for the spatial and grain size 
distributions. It accounts for the spatial filtering of the in- 
terferometers and uses the Mie theory to self-consistently 
simulate anisotropic scattering depending on grain proper- 
ties. Considering a wide range for the dust parameters, we 
produced disc models that fit well our dispersed H-band 
best-fit flux ratios, the K-band contrast computed from 
archival VLTI/VINCI measurements, and upper limits in 
the mid-infrared provided bv single-dish telescopes f Telesco 
et al. I20050 . The results show that light scattered in the 
known debris disc might account for a significant part (up 
to 70%) of the near-infrared excess flux, while producing a 
very good match to archival mid- to far-IR data through 
its thermal emission. This maximum value of 70% must 
be considered as very conservative since the Mie theory 
tends to overestimate forward- and back-scattering owing 
to the spherical grain hypothesis. The spectral slope in the 
near-IR is reproduced well by this scenario, but additional 
hot circumstellar material must be present in the vicinity of 
l3 Pic in order to account for the remainder of the measured 
H-band flux ratio. 

The presence of hot material around (3 Pic would not 
be a surprise and might have various origins, such as stel- 
lar wind (free- free emission), mass- loss events, hot gas, or 
hot exozodiacal dust. Stellar winds and mass- loss events 
can, however, be excluded with very good confidence since 



they are expected to be very weak in t he case of A-type 
stars (see discussion in lAbsil et al.|[2008h . Hot gas is likely 
to b e present in the inner system (< a few AU. Hobbs 
et al. Il988^ and must emit in the continuum in order to 
account for our measurements. Assuming that the gas is 
produced by evaporating comets and brak ed by a ring of 
neutral hydrogen located at about 0.5 AU. 'Lagrange et al.l 
(1998} derived temperature estimates of 1500-2000 K with a 
gas production rate of ~10~^^ Mq/yt. The question is thus 
whether hot gas could accumulate in sufficient quantities to 
contribute (at least partially) to the detected near-infrared 
emission. This is, however, beyond the scope of this paper 
and will not be addressed here. A more likely scenario is the 
presence of hot exozodiacal dust as proposed for ol der A- 
tvpe stars fe.g.. lAbsil et"aIll2006l ; [Defrere et al.ll2011l . Men- 
nesson et al., submitted to ApJ). In the case of /? Pic, the 
hot dust in the inner region of the disc could be produced 
by comets ("falling evaporating bodies'') dynamica lly per- 
turbed by one or several planets (jBeust et al.lll990l . Bonsor 
et al., submitted to A&A). More investigations are clearly 
needed to establish the most realistic scenario and, particu- 
larly, the balance between hot material and scattered light. 
This will be the subject of a forthcoming paper (Lebreton 
et al., in prep) that will include both near- infrared (H and 
K bands, this work) and mid-infrared (di Folco et al., in 
prep.) resolved observations. 
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